Spin-spin correlators in Majorana representation 
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In the Majorana representation of a spin 1/2 we find an identity which relates spin-spin correlators 
to one-particle fermionic correlators. This should be contrasted with the straightforward approach 
in which two-particle (four-fermion) correlators need to be calculated. We discuss applications to the 
analysis of the dynamics of a spin coupled to a dissipative environment and of a quantum detector 
performing a continuous measurement of a qubit's state. 
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I. INTRODUCTION 

An analysis of spin dynamics involves calculations of 
spin correlators. Spin operators do not satisfy the Wick 
theorem, and various methods have been used to still 
enable the use of perturbative (diagrammatic) methods. 
One of the approaches is based on the Majorana-fermion 
representation of spin operators. This approach has 
a long history and was applied recently to condensed- 
matter problems 

In this approach one introduces three Majorana 
fermions r\ x ,y,z (per spin) and expresses the spin (or 
Pauli) operators via these fermions: 



Another Dirac fermion g = c z is still needed to construct 



Oy = -ir) z r) x 



(1) 



Obviously, an analysis of spin-spin correlations based on 
Eq. Q requires calculations of a four-fermion correlator. 
In this note we show that spin correlators coincide with 
certain two-fermion (i.e., one-particle) correlators. We 
discuss under which conditions this identity may simplify 
the analysis and discuss its particular applications. 



II. MAJORANA REPRESENTATION 

Majorana fermions satisfy the anti-commutation rela- 
tions: r] a r]{3 — —rjprja for a ^ (3 and rfe = 1, and are 



real: ry a = r\ a (a, (3, 7 = x, y, z). These properties ensure 
that the representation reproduces the commutation 
relations of the spin algebra. An important feature of the 
representation is the fact that the spin operators are 
bilinear ('bosonic') combinations of Dirac (annihilation / 
creation) fermionic operators. 

One may construct the three Majorana fermion oper- 
ators (for a spin) out of three different Dirac fermions, 
Va = c a + c a , each annihilation operator acting in its 
own two-dimensional Hilbert space. This ensures the an- 
ticommutation relations and the property rfc = 1. The 
whole Hilbert space is then 8-dimensional. 

The dimensionality of the Hilbert space may be re- 
duced if two Majorana fermions, e.g. rj x and r] y , are con- 
structed as linear combinations of a single Dirac fermion 
/ and its conjugate r] x = / + P and r] y = i(p - /). 



the third Majorana fermion r\ z = g 
Majorana-Dirac picture 

0+ = Vzf 
v- = PVz 
<7, = 1 - 2/V 
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In this mixed 



(2) 



This 'drone-fermion' representation was used for the 
analysis of magnetic systems in the 60s (see RefsA^Sii). 
The whole Hilbert space in this representation is 4- 
dimensional. Depending on the rotational symmetry, this 
representation may be more convenient than 

We have described this construction in the Hilbert 
space of two Dirac fermions / and g. Alternatively one 
can view it as two replicas of the original spin. Indeed, 
let us label the basis states in the following way: de- 
note the state without /- and g-fermions by |f ) = 1 00) , 
and also |t 6 > = |01) = ff t | 0), || 6 ) = |10) = /t |00), 
and \i a ) = |11) = /V |00>. For the state \s n ) the first 
index s =f / J, denotes the spin component, while the 
second n = a/b labels the spin copy. One notices that 
the spin operators c+ , er_ , and a z do not mix the a- and 
6-subspaces, i.e., they operate in the same way on the 
two 'copies' of the spin. Further, one may view the index 
n = a/b as an isospin and introduce the respective Pauli 
isospin operators t x . v ,z- I n particular, t x is the 'copy- 
switching' operator: t x \s a ) = \s b ) and t x \s b ) = \s a ). 
The fermionic operators can be expressed as 



/ 



>+ix, 



(3) 



Further, r] a = cr a T x , for any a = x, y, z. Accordingly, 
in terms of the fermionic operators we have t x = (1 — 
2/ /) Vz = ^iilxVyVz- The operator t x commutes with 
'all other' operators: with <j + , <t_, <j z , /, and r\ z . 



III. REDUCTION OF THE SPIN-SPIN 
CORRELATORS 

A physical Hamiltonian depends on the spin operators 
a x , o yi a z (and on other degrees of freedom, e.g., the 
electrons in the Kondo problem). Thus the operator t x 
commutes with any Hamiltonian and we obtain 



= (T x (t)r) a (t)T x (t')r)p(t')) 



iv a (t) Vp (t'))) . (4) 
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Here we have used the fact that t x is time- independent, 
commutes with the Majorana fcrmions, and that t£ = 1. 
For example, we obtain 

(o- x (t)a x (t')) = {r) x {t)r]x{t')) 

= ([fit) +fHm w) +/ + (*')]> • (5) 

In certain situations the identity may simplify the 
calculations. Indeed, if the Majorana representation Q, 
(0) is used for a calculation of spin-spin correlations, a 
four-fermion correlator needs to be evaluated. In a typi- 
cal situation, the lowest-order contribution is given by a 
loop-like diagram which involves two fermionic propaga- 
tors. A straightforward evaluation of higher-order contri- 
butions requires an analysis of the self-energy corrections 
to these propagators as well as of the vertex corrections. 
The use of the relation Q reduces the task to the eval- 
uation of a single self-energy. 

In general, evaluation of this self-energy to all orders 
in the perturbative expansion involves complicated dia- 
grams and, in particular, other self-energies and vertex 
corrections. Hence an involved calculation may still be 
needed. Nevertheless, the relation (@J may be useful if 
the needed self-energy part can be estimated reliably, for 
instance, by calculating the low-order contributions. We 
discuss two examples in Section Ivl 

Note also that the relation Q may be straightfor- 
wardly generalized to multi-spin correlators, the evalu- 
ation of which then reduces to multi-fermion correlators. 
Further, in a problem, that involves many spins repre- 
sented via Majorana fermions, (e.g., a lattice spin model) 
the r^-operators for different spins (sites) may be in- 
volved in a relation similar to Q . In this case their prod- 
uct does not drop out of the calculation (unlike = 1 
for one spin). 

IV. GAUGE-INVARIANCE CONSIDERATIONS 

One might be concerned by the fact that Eq. (J3J re- 
duces a correlator of physical quantities to a correlator of 
"unphysical" operators. Another way of expressing this 
concern is to invoke the gauge symmetry. Indeed the Ma- 
jorana representations and (J2J possess the discrete Z?- 
symmetry rj a — > —rj a . As an operator which realizes the 
symmetry transformation one may choose r y — i(g* — g) 
(cf. Ref£). For example, r\ z — » T y r\ z Ty X = —rj z . Thus, 
the forth Majorana fermion allowed in the Hilbert space 
of the representation J2J) generates the symmetry trans- 
formation. One can as well use as a generator r z , which 
just flips the sign of the wave function of the 6-spin, keep- 
ing that of the a-spin intact. 

Consider now a time-dependent gauge transformation, 
which transforms a wave function |^) to U \ where 

(Z7T \ 
y t v W)) (6) 

and <f>(t) =0,1. It transforms the operators r/ Q — > 
(-l) 0(t) ?7 Q , t x -> (-l)^*)r x . Thus the operator t x is 



now time-dependent and no longer commutes with the 
Hamiltonian. Indeed, the gauge-transformed Hamilto- 
nian reads 

H = UHlf- 1 +iUU- 1 = H -^<p{t) Ty , (7) 

and the last (gauge) term does not commute with t x . 
Thus we find that 

(Va(t)Va(t')) - (-lr^-^'HVa^Vait')) , (8) 

i.e., the single- fermion correlators are not gauge- 
invariant. 

The relation J3J is written in a fixed gauge, in which 
4>{t) = 0. In an arbitrary gauge, the factor (— l)^( t ) _< ^(* ) 
should be added on the rhs making the identity gauge- 
invariant. In path- integral calculations the integration 
is performed over all configurations (all gauges). If a 
saddle-point solution breaks the .^-symmetry (and hence 
has a counterpart) fluctuations which involve switching 
between two mean-field solutions play an important part 
and need special attention. 8 However, for perturbative 
calculations around a trivial saddle point, e.g. above the 
Kondo temperature, the 'fixed-gauge' relation Q might 
be quite useful. 

V. APPLICATIONS 

The relation J3J can simplify calculations, since to eval- 
uate a single-fermion Green function one needs only to 
evaluate a self-energy. In a perturbative regime, when 
one may restrict oneself to evaluating the lowest-order 
contribution, Eq. J3J simplifies the task. 

One example of such a problem is dicussed in Ref. 9 
where continuous measurement of the state of a spin 
(qubit) by a quantum detector is analyzed and the ouput 
noise of the detector as well as the spin dynamics is 
studied. A calculation of a spin-spin correlator in the 
Majorana representation requires an evaluation of a two- 
fermion loop and involves the analysis of two self-energy 
parts and a vertex correction. It turns out that the ver- 
tex correction is important already in the lowest order 
and may even cancel one of the self-energy contributions. 
This fact may be understood and the calculations are 
considerably simplified if the relation (JIJ is invoked. The 
obtained result coincides with an alternative calculation, 
in which instead of using the Majorana representation 
to enable the use of the Wick theorem one follows the 
dynamics of the spin directly. This approach-^ may be 
useful if the problem involves only a small number of 
spins (or other 'non-Wick' degrees of freedom). 

Here we consider another example: the dissipative dy- 
namics of a spin coupled to an environment, for instance, 
the spin-boson model. Consider the case of a purely 
transverse coupling: 

H = --B <j z - ±X(t) <j x + tfbath , (9) 
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where X is a fluctuating bosonic observable of the bath, 
whose Hamiltonian -ffbath determines the statistics of 
fluctuations. We consider gaussian, but not necessarily 
equilibrium, fluctuations. Using the Keldysh technique 
we perform a calculation considering the bath-spin cou- 
pling (the second term on the rhs of Eq. 0) as a pertur- 
bation. 

Below we use the notations of Refiii. Let us introduce 
the matrix Green function of the bath: 



• • ■ • 



E = 



G x = -i(T K X(t)X(t')) 



Its Keldysh component 



f~<K /"-r> 



G< = -2iS x 



(10) 



(11) 



and the difference between the retarded and advanced 
components 



V 



FIG. 1: The lowest order non-vanishing contributions tor the 
self-energies and ~S V - The solid lines denote the propaga- 
tors of the bath (X), the dashed line — of the /-fermion, and 
the dotted lines — of the Majorana fermion r\. (See text.) 



G 



x 



Gx 



G x 



-2iA 



x 



(12) 



are related to the symmetrized and antisymmetrized 
correlators: (X(t)X(0)) = S x (t) + A x (t), S x (-t) = 
Sx(t), Ax{— t) = — A x (t), In equilibrium they are re- 
lated by the fluctuation-dissipation theorem: Sx{^>) — 
Ax(uj)coth{huj/2k B T). Here G x = -i(X{t)X(0)) and 
G<= r i{X(0)X{t)). 

Similarly, for the Majorana fermion ?; = rj z we define 
G v = —i(TnT]{t)ri(t')) . The bare Green functions are 
G^q = —i and G< = i. For the /-fermion we use the 
Bogolubov-Nambu spinors * = (/, p) T and '3> t = (/ T , /) 
and define a matrix G* = -i{T K ^(t)^ (*'))• 

Calculation of the spin propagators (in a stationary 
state) reduces due to the relation (0} to the evaluation of 
the fermion Green functions G^, and G^. These functions 
can be found from the the Dyson (kinetic) equations: 



G,i 



G 



*,o 



G-q — G, h o ^"n , 



(13) 



(14) 



where E^ and are the self-energies. All the quan- 
tities in Eq. (| 1 31) are 4x4 matrices (with the Nambu 
and Keldysh components). In the frequency domain the 
operator G^o is given by 



w-B 

i i uj + B 

*.° I lu — B 

uj + B 



(15) 



The lowest-order contributions to the self-energies 

^< = I ~7 v'i ) ■ = I Z'\ I (i7) 



o 



yA 



E^ 



are shown in Fig. ^ (the matrix structure after the 
Keldysh rotation is giveniA*i£). We find that = 



(i/4)AG>G> and E< = (i/4)AG<G<, where A = 



?> r. 
1 -1 



-1 1 



ia a matrix in the Nambu space. Sim- 



ilarly, we obtain E> = (i/4) G x ( 1 -1 ) G> ( _j 



and E< = 



(i/4)G< (1 -1)G< 



While the 



retarded and advanced components of the bare fermion 
Green functions are known, its Keldysh component con- 
tains information about the distribution function and is 
sensitive to the environment (cf. Refii^). 

To perform the calculation, one introduces 
the distribution function h v via G^(w) = 

h v (w) (G*f(w) - G£(uj)). Evaluating E* from the 
first diagram in Fig. ^ we find 

E£(w)-E£( W ) - ^ [G x (u>) + h n (0)(G x (u) G x (cj))} , 

(18) 

s f M = ^A [(G x (u) - G x {uj)) + hr,(Q)Gx(u)] . (19) 

The symmetry relation h v (—uj) = —h v (uj) implies 
h v (0) = 0. Hence without evaluating we find 



while for G^.q we obtain 



r -i_(u}/2 



(16) 



We disregard infinitesimal imaginary terms in Eqs. (|15l 
11611 since they are superseded by the self-energy parts in 
Eqs. JHEJl. 



and 



Sf ( w ) - E^M 



•-ASx(w) 



\A x {io) 



(20) 



(21) 



The real parts of the retarded and advanced self-energies 
give the non-equilibrium generalization of the Lamb shift, 
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i.e., renormalize the level splitting B. This renormaliza- 
tion is small if the noise is weak and has a non-singular 
spectrum. 

Substituting the self-energy into the Dyson equa- 
tion 113(1 we obtain 

G-i (~*~ i 

/ -ir(u) iT(oj) -iA x (w)/2 iA x (w)/2 

iT(u>) -ir(w) iA x {w)/2 -iA x {w)/2 

"00 iT(uj) -iF(u) 

\ -iT(u>) iF(u) 



(22) 



where T(ui) = S x (u})/4. Inverting Eq. 1)22(1 we obtain 



n R/A _ 



uj + B± iT(uj) ±iT{uj) 
±iT(ui) uj-B±iT(oj) 

(to 2 - B 2 )±2iuT(u;) 



iA x {») ( -(u + B) 2 lu 2 -B 2 
{ uj 2 -B 2 -{uj~B) 2 

(to 2 - B 2 ) 2 + 4w 2 i» 2 



(23) 



(24) 



From Eqs. ((23(1 and 124|) we see that, at least in this 
order of perturbation theory, we can use the relation 
Gf(w) = h 9 (u)(G^(u)-Gi(u)), where h*(-w) = 
— h^(uj) and (for any pair ij of Nambu indices) 



= A x {lo) 



(25) 



a quantity of the zeroth order in the coupling constant. 
Using these results we can calculate £,,: 

E«M-£>) = i[G$( W -i?) 

+/i*(B)(G£(o; - S) - G x {lo - B))] 
+[B^-B], (26) 

E£( W ) = i[(G£( W -fl)-G£( W -B)) 

+ h< s ,(B)G§ {uj - B)) + [B -> -B] . (27) 

We evaluate the self-energy E,, near the pole ui = of 
the Green function G„: 



where the bosonic time-ordering is chosen for the spin 
components. From Eq. © we deduce that 



XX 



Since G^ 



-i(a x (t)a x (t')) = ( 1 1)G>(}) . (31) 



2 t^* 



G^], we obtain 



(SxM + AxM)B 2 
(w 2 -B 2 ) 2 +4cj 2 r(w) 2 



(32) 



This result describes the peaks in the spectrum of spin 
corelations at us = ±B, with the width corresponding 
to the dephasing rate T^ 1 = T(B) = S x (B)/4, thus 
reproducing the form known, e.g., from the analysis of 
the Bloch-Rcdfield equations 

Similarly, for the Green function of the Majorana 
fermion r\ we obtain 



where 



t _ Sx(B) 



uj±2iT 



A\{B) 
S\{B) 



The Keldysh component vanishes, G^ 



(33) 
(34) 

0, and we find 



»n> = (a z (t)a z (t')) u 



iG> 



4r 



4r 2 



(35) 



This may be compared to the peak shape obtained in 
the Bloch-Redfield approach^ It can be expressed in 
terms of the average value (a z ) = A X (B)/S X (B) of the 
2-component: 



(a z ) 2 2ttS(u;) + [1 - (a z ) 2 



217 1 
J 2 + Tf 2 



(36) 



where the relaxation time is given by T^ 1 — 2T 2 -1 . 
The result <(33)l reproduces the high-frequency behav- 
ior of Eq. 1(36(1 and the 'width of the peak', which may 
be read off from the high-w asymptotics, iH zz ~ 2(1 — 
(cz) 2 )/(TiLli 2 ). This indicates that at lower frequencies 
higher-order corrections (contributions to £,,) are impor- 
tant, which requires a further analysis. 



Ef(0)=0 



and 



£*(0)-S» = -~S x (B) 1- 



A\{B) 
S\{B) 



(28) 



(29) 



Now, with the acquired knowledge of single-fermion 
Green functions, we can evaluate various spin-spin cor- 
relators. We start with 



n. 



-i(T K a x (t)a x (t'))) , 



(30) 



VI. SUMMARY 

In summary, we discussed an identity relating spin 
correlations to single-particle Majorana-fermion propa- 
gators and problems where the use of this relation sim- 
plifies the analysis. 
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tions. Y.M. was supported by the Humboldt foundation, 



5 



the BMBF, and the ZIP Programme of the German gov- of this manuscript a preprint ^ appeared which partially 
crnment. overlaps with our work. 

Note added: At the final stage of the preparation 



1 A. M. Tsvelik, Phys. Rev. Lett. 69, 2142 (1992). 

2 P. Coleman, E. Miranda, and A. M. Tsvelik, Phys. Rev. 
Lett. 70, 2960 (1993). 

3 A. M. Tsvelik, Quantum Field Theory in Condensed Mat- 
ter Physics (Cambridge Univ. Press, Cambridge, 1995). 

4 D. C. Mattis, Theory of Magnetism (Harper and Row, New 
York, 1965). 

5 R. P. Kenan, J. Appl. Phys. 37, 1453 (1966). 

6 H. J. Spencer and S. Doniach, Phys. Rev. Lett. 18, 994 
(1967). 

7 H. J. Spencer, Phys. Rev. 171, 515 (1968). 

8 P. Coleman, L. B. Ioffe, and A. M. Tsvelik, Phys. Rev. B 
52, 6611 (1995). 



9 A. Shnirman, D. Mozyrsky, and I. Martin, in preparation 
(2003). 

H. Schoeller and G. Schon, Phys. Rev. B 50, 18436 (1994). 

1 J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

2 A. Kamenev, cond-mat/0109316 (2001). 

3 O. Parcollet and C. Hooley, Phys. Rev. B 66, 085315 
(2002). 

4 See, e.g., A. Shnirman, Yu. Makhlin, and G. Schon, Phys. 
Scr. T102, 147 (2002) and references therein. 

5 W. Mao, P. Coleman, C. Hooley, and D. Langreth, cond- 
mat/0305001 (2003). 



